Strange Nonchaotic Bursting in the Quasiperiodically Forced Hindmarsh-Rose Neuron 
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We study the transition from a silent state to a bursting state by varying the dc stimulus in the 
Hindmarsh-Rose neuron under quasiperiodic stimulation. For this quasiperiodically forced case, a 
new type of strange nonchaotic (SN) bursting state is found to occur between the silent state and 
the chaotic bursting state. This is in contrast to the periodically forced case where the silent state 
transforms directly to a chaotic bursting state. Using a rational approximation to the quasiperiodic 
forcing, the mechanism for the appearance of such an SN bursting state is investigated. Thus, a 
smooth torus (corresponding to a silent state) is found to transform to an SN bursting attractor 
through a phase-dependent subcritical period-doubling bifurcation. These SN bursting states, to- 
gether with chaotic bursting states, are characterized in terms of the interburst interval, the bursting 
length, and the number of spikes in each burst. Both bursting states are found to be aperiodic com- 
plex ones. Consequently, aperiodic complex burstings may result from two dynamically different 
states with strange geometry (one is chaotic and the other one is nonchaotic). Thus, in addition to 
chaotic burstings, SN burstings may become a dynamical origin for complex physiological rhythms 
which are ubiquitous in organisms. 

PACS numbers: 05.45. Ac, 05.45.Df, 87.19.L- 



I. INTRODUCTION 

To probe dynamical properties of a system, one often 
applies an external stimulus to the system and study its 
response. Particularly, periodic stimulation to biological 
oscillators has attracted much attention in various sys- 
tems such as the embryonic chick heart-cell aggregates 
[H, the squid ^nt axon 0, and the cortical pyra- 
midal neurons [4]. Rich regular (such as phase locking 
and quasiperiodicity) and chaotic responses were found 
in these periodically forced systems [1, @]. In contrast, 
quasiperiodically forced case has received little attention 
[Ti, 8], and hence further intensive investigation on dy- 
namical responses of quasiperiodically forced biological 
oscillators is necessary. 

Here, we are interested in neural bursting activity [al- 
ternating between a silent phase and an active (burst- 
ing) phase of repetitive spikings] [ol]. Cortical intrin- 
sically bursting neurons, thalamocortical relay neurons, 
thalamic reticular neurons, and hippocampal pyramidal 
neurons are representative examples of bursting neurons 
[l0| . We are particularly concerned about dynamical 
responses of bursting neurons subject to quasiperiodic 
stimulation. Strange nonchaotic (SN) attractors typi- 
cally appear in quasiperiodically forced dynamical sys- 
tems |ll| - [l4| . They exhibit some properties of regular as 
well as chaotic attractors. Like regular attractors, their 
dynamics is nonchaotic in the sense that they do not 
have a positive Lyapunov exponent; like usual chaotic 
attractors, they have a geometrically strange (fractal) 
structure. Hence, SN burstings are expected to occur 
in quasiperiodically forced bursting neurons. 
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This paper is organized as follows. In Sec. |TI1 we con- 
sider the Hindmarsh-Rose (HR) neuron model for burst- 
ing neurons which was originally introduced to describe 
the time evolution of the membrane potential for the 
pond snail and investigate the transition from 

a silent state to a bursting state by varying the dc stim- 
ulus. This work is in contrast to previous works on the 
effect of the quasiperiodic forcing on the self-oscillating 
neurons in the spiking state of self-sustained oscillations 
of the membrane potential Q . In the periodically forced 
case {i.e., in the presence of only one ac stimulus source), 
an intermittent transition from a silent state (with sub- 
threshold oscillations) to a chaotic bursting state occurs 
when the dc stimulus passes a threshold value. Effect 
of the quasiperiodic forcing on this intermittent route 
to chaotic bursting is investigated by adding another in- 
dependent ac stimulus source. Thus, unlike the case of 
periodic stimulus, a new type of SN burstings are found 
to occur between the silent state and chaotic bursting 
state as intermediate ones. Using a rational approxima- 
tion to the quasiperiodic forcing (TJ, , we investigate 
the mechanism for the appearance of such SN burstings. 
Thus, a smooth torus, corresponding to a silent state, 
is found to transform to an SN bursting attractor via 
a phase-dependent subcritical period-doubling bifurca- 
tion. Together with chaotic burstings, these SN burstings 
are characterized in terms of the interburst interval, the 
bursting length, and the number of spikes in each burst. 
Both the chaotic and SN bursting states are found to 
be aperiodic complex ones. Such aperiodic complexity 
comes from the strange geometry of both bursting states 
with qualitatively different dynamics (one is chaotic and 
the other one is nonchaotic) . We note that complex phys- 
iological rhythms, which are central to life, are ubiquitous 
in organisms 6|. Hence, in addition to chaotic burstings, 
SN burstings may also serve as a dynamical origin of such 
complex bodily rhythms. Finally, a summary is given in 
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II. SN BURSTINGS IN THE 
QUASIPERIODICALLY FORCED HR NEURON 

We consider a representative HR bursting neuron 
model @, which is quasiperiodically forced at two 
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where /cxt = -^dc + sin(27r/it) + A2 sin(27r/2i), a = 1, 
b = 3, c — I, d = 5, s — 1, r — 0.001, and xq = —1.6. 
Here, t is the time [measured in units of miUisecond (ms)], 
X is the membrane potential variable, y is the recovery 
variable, z is the slow adaptation current, Idc is a dc stim- 
ulus, Ai and A2 are amplitudes of quasiperiodic forcing, 
and u){= /2//1) is irrational (/i and f2- measured in 
units of kHz) . 

To obtain the Poincare map of Eq. ([1]), we make a 
normalization fit — > t, and then Eq. ^ can be reduced 
to the following differential equations: 
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where x = (z, y, z) and /oxt = ^dc + Ai sin(27ri) + 
sin(27r6'). The phase space of the quasiperiodically 
forced HR oscillator is five dimensional with coordinates 
x, 2/, z, 0, and t. Since the system is periodic in and t, 
they are circular coordinates in the phase space. Then, 
we consider the surface of section, the x-y-z-9 hypersur- 
face dX t — n (n: integer). The phase-space trajectory 
intersects the surface of section in a sequence of points. 
This sequence of points corresponds to a mapping on the 
four-dimensional hypersurface. The map can be com- 
puted by stroboscopically sampling the orbit points v„ 
[= (x„. On)] at the discrete time n (corresponding to mul- 
tiples of the first external driving period Ti). We call the 
transformation v„ — > v„_(_i the Poincare map, and write 
v„+i = P(v„). 

Numerical integration of Eqs. ^ and © is done us- 
ing the fourth-order Runge-Kutta method. Dynamical 
analysis is performed in both the continuous-time system 
(ie., flow) and the discrete-time system (i.e., Poincare 
map). For example, the time series of the membrane po- 
tential x{t), the phase flow, the interburst interval, the 
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FIG. 1: Intermittent transition to a chaotic bursting state 
for the case of periodic forcing with Ai = 0.5 and /i = 30 
Hz {A2 = 0). (a) Time series of x(t) and (b) projection of 
the phase flow onto the x — z plane for the silent state ex- 
hibiting subthreshold oscillations for /dc = 0.3. (c) Time 
series of x{t) and (d) projection of the phase flow onto the 
X — z plane for the chaotic bursting state for /dc = 0.5. 
(e) Bifurcation diagram {i.e., plot of x versus /dc) in the 
Poincare map. We obtain attractors by iterating the Poincare 
map at the 500 equally-spaced values of /dc in the range of 
/dc £ [0.2, 0.57]. For each chosen /dc, we choose a random ini- 
tial point (a;(0), 2/(0), z(0), 0(0)) with uniform probability in 
the range of a;(0) G (-2,2), y(0) G (-16,0), 2(0) £ (0,0.4), 
and 9{Q) G [0, 1), and obtain the attractor through the 200- 
times iterations of the Poincare map after the transients of 
the 1000 Poincare maps. 



bursting length, and the average number of spikes per 
burst are obtained in the flow. On the other hand, the 
Lyapunov exponent [17] and the phase sensitivity expo- 
nent [1^ of an attractor are calculated in the Poincare 
map. To obtain the Lyapunov exponent of an attrac- 
tor in the Poincare map, we choose 20 random initial 
points {{xi{Q), j/,(0), Zi(0), 6*^(0)); i = 1, . . . , 20} with uni- 
form probability in the range of Xi{Q) G (—2,2), yi{Q) G 
(-16,0), z,(0) G (0,0.4), and 6',(0) G [0, 1)- For each ini- 
tial point, we get the Lyapunov exponent [l7[, and choose 
the average value of the 20 Lyapunov exponents. (The 
method of obtaining the phase sensitivity exponent will 
be explained below.) 

Here, we set u to be the reciprocal of the golden mean 
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[i.e., (jj = (-s/S — l)/2], and numerically investigate dy- 
namical transition from a silent state to a bursting state 
by varying Idc in the HR neuron under external stimu- 
lus. We first consider the case of periodic forcing {i.e., 
A2 = 0) for Ai = 0.5 and /i = 30 Hz. Figures [^a) 
and[IJb) show the time series of x{t) and the projection 
of the phase flow onto the x-z plane for the silent state 
when /dc = 0.3. We note that this silent state with the 
largest Lyapunov exponent ci ~ —0.133 exhibits sub- 
threshold oscillations. As /dc passes a threshold value of 
-^dc — 0.416 721, a chaotic bursting state appears. Burst- 
ing activity [alternating between a silent phase and an 
active (bursting) phase of repetitive spikings] of the HR 
neuron is shown in Fig. [Ijc) for /dc — 0.5. This kind 
of bursting occurs on a chaotic hedgehoglike attractor 
with (Ti ~ 0.406 [the body (spines) of the hedgehoglike 
attractor corresponds to the silent (bursting) phase], as 
shown in Fig.[T]Jd). This transition from a silent state to a 
chaotic bursting state is investigated by varying /dc in the 
Poincare map. Figure[TJe) shows the bifurcation diagram 
(i.e., plot of X versus /dc)- The solid curve represents a 
stable fixed point corresponding to the silent state. The 
stable fixed point loses its stability for /dc = 0.416 721 
via a subcritical Hopf bifurcation when a pair of complex 
conjugate stability multipliers passes the unit circle in 
the complex plane, and then a chaotic bursting attrac- 
tor, corresponding to a chaotic bursting state, appears. 

From now on, we consider the quasiperiodically forced 
case for Ai = 0.5 and /i = 30 Hz. Each state is char- 
acterized by both the largest (nontrivial) Lyapunov ex- 
ponent (Ti, associated with dynamics of the variable x 
[besides the (trivial) zero exponent, related to the phase 
variable 9 of the quasiperiodic forcing] and the phase sen- 
sitivity exponent 5. The exponent 5 measures the sensi- 
tivity of the variable x with respect to the phase of the 
quasiperiodic forcing and characterizes the strangeness of 
an attractor [31 . A (regular) silent state has a negative 
Lyapunov exponent (i.e., cti < 0) and has no phase sensi- 
tivity (i.e., 5 — 0). On the other hand, a chaotic bursting 
state has a positive Lyapunov exponent cti > 0. In ad- 
dition to them, a new type of SN bursting states that 
have negative Lyapunov exponents (cti < 0) and positive 
phase sensitivity exponents {5 > 0) appear. Due to their 
high phase sensitivity, SN bursting states have a strange 
fractal phase space structure. For small A2, a direct tran- 
sition from a silent state to a chaotic bursting state oc- 
curs, as in the periodically forced case of A2 = 0. As an 
example, we consider the case oi A2 = 0.2 where a transi- 
tion to chaotic bursting occurs for /dc = 0.3963. Figures 
[2)Ja) and [IJb) show a smooth torus with cri ~ —0.036 
(corresponding to a silent state) and a chaotic bursting 
attractor with cti ~ 0.154 (corresponding to a chaotic 
bursting state) for /dc — 0.39 and 0.4, respectively. Fig- 
ure lUc) shows the Lyapunov-exponent diagram (i.e., 
plot of cTi vs. Idc). As /dc is increased to the transition 
point, (Ti of the smooth torus increases to zero, and then 
a chaotic bursting attractor with a (finite) positive cti ap- 
pears (i.e., a finite jump for the value of <ti seems to oc- 
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FIG. 2: Transition from a silent state to a bursting state for 
the quasiperiodically forced case when Ai = 0.5 and /i = 30 
Hz. We set A2 = 0.2 in (a)-(c). Projections of attractors 
onto the 9 — x plane in the Poincare map are shown for the 
(a) silent and (b) chaotic bursting states when /dc ~ 0.39, 
and 0.4, respectively, (c) Lyapunov-exponent diagram (i.e., 
plot of o"! versus Idc). For each randomly chosen initial point, 
we get the Lyapunov exponent by following the 10* Poincare 
maps after the transients of the 10* Poincare maps. We set 
A2 — 0.5 for (d)-(i). Time series of x{t) of (d) the silent state 
for /dc = 0.21, (e) the SN bursting state for /dc = 0.24, and 
(f) the chaotic bursting state for Idc ~ 0.29. Projections of 
the phase flows onto the x — z plane for (g) the silent state 
for Idc = 0.21, (h) the SN bursting state for Idc = 0.24, and 
(i) the chaotic bursting state for Idc = 0.29. 



cur). However, for A2 larger than a threshold A2 (~ 0.4), 
SN bursting states appear between the silent and chaotic 
bursting states. As an example, we consider the case 
of A2 = 0.5 and investigate dynamical behaviors of the 
quasiperiodically forced HR neuron by varying /dc. As 
/dc passes a threshold I^^ (~ 0.2236), the silent state be- 
comes unstable, and a transition to an SN bursting state 
occurs. As /dc is further increased and passes another 
threshold value of /dc — 0.2703, the SN bursting state 
transforms to a chaotic bursting state. Figures [2l|d)-[2Uf) 
show the time series of the membrane potential x{t) of 
a silent state (exhibiting subthreshold oscillations), an 
SN bursting state, and a chaotic bursting state for /dc — 
0.21, 0.24, and 0.29, respectively. For these three cases, 
the phase flows of the silent, SN bursting and chaotic 
bursting states are also given in Figs. [Hg)-llUi), respec- 
tively. 

The silent and bursting states for A2 — 0.5 are ana- 
lyzed in terms of the largest Lyapunov exponent cti and 
the phase sensitivity exponent S in the Poincare map. 
Projections of attractors onto the — x plane for /dc = 
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FIG. 3: Characterization of silent and bursting states for 
A\ = 0.5, A-i, — 0.5, and and /i = 30 Hz. Projections 
of attractors onto the 9 — x plane in the Poincare map are 
shown for the (a) silent, (b) SN bursting, and (c) chaotic 
bursting states when Idc ~ 0.21, 0.24, and 0.29, respectively, 
(d) Lyapunov-exponent diagram (i.e., plot of ai versus /dc); 
(Ti for the SN bursting attractor is shown in black. For each 
randomly chosen initial point, we get the Lyapunov exponent 
by following the 10* Poincare maps after the transients of the 
10^ Poincare maps, (e) Phase sensitivity functions pj^' are 
shown for the silent and SN bursting attractors when Jdc=0.21 
and 0.24, respectively. For the case of SN bursting attractor, 
the graph is well fitted with a dashed straight line with slope 
S ~ 1.74. (f) Plot of the phase sensitivity exponent S versus 
A/dc (= ^dc — -^dc) for the SN bursting attractor; I^^ ~ 0.2236. 



0.21, 0.24, and 0.29 are given in Figs. [2Ia)-[3Uc), respec- 
tively. For the silent smooth torus exists in the 
6 — x plane [see Fig.[3][a)]. On the other hand, nonsmooth 
bursting attractors appear for both SN and chaotic burst- 
ing states, as shown in Figs. [3l^b) and|31[c). A dynam- 
ical property of each state is characterized in terms of 
the largest Lyapunov exponent cri (measuring the de- 
gree of sensitivity to initial conditions) . The Lyapunov- 
exponent diagram is given in Fig. El^d) . When passing the 
bursting transition point I^^i— 0.2236), an SN bursting 
attractor appears. The graph of ai for the SN bursting 
state is shown in black, it is nearly flat, and its value 
is negative as in the case of smooth torus. However, as 
Idc passes the chaotic transition point I^c (— 0.2703), 
a chaotic bursting attractor with a positive ai appears. 



Unlike the case of direct transition from a smooth torus to 
a chaotic bursting attractor, tJi seems to increase contin- 
uously from zero without jump [compare Fig. ^d) with 
Fig. d^c)]. Although SN and chaotic bursting attractors 
are dynamically different, they both have strange geom- 
etry leading to aperiodic complex burstings. To charac- 
terize the strangeness of an attractor, we investigate the 
sensitivity of the attractor with respect to the phase 9 of 
the external quasiperiodic forcing |13| . This phase sen- 
sitivity may be characterized by differentiating x with 
respect to 6* at a discrete time t = n. Using Eq. 
we may obtain the following governing equation for 
(z = 1,2,3), 



d / dxi \ \ - dFi dxj dFi 



.1 



80 80 



(3) 



where (xi, 2:2, 2^3) = {x,y,z) and Fi's {i = 1,2,3) 
are given in Eq. Starting from an initial point 

(x(0), ^(0)) and an initial value dx/80 = for t = 0, we 

may obtain the derivative values of Sn^ (= dxi/d0) at 
all subsequent discrete time t = nhy integrating Eqs. ^ 

and ([3|). One can easily see the boundedness of Sn^ by 
looking only at the maximum 



7^^(x(O),0(O)) 

(4) 

We note that 7^^(x(0), 0(0)) depends on a particular tra- 
jectory. To obtain a "representative" quantity that is in- 
dependent of a particular trajectory, we consider an en- 
semble of randomly chosen initial points {(x(0), 6'(0))}, 
and take the minimum value of with respect to the 
initial orbit points 13], 



max |5«(x(O),0(O))| (z = 1,2,3). 

0<n<A' 



■ N 



min i^\x{Q),0{Q)) [i 
{(x(o),e(o))} 



1,2,3). (5) 



Figure[3fe) shows a phase sensitivity function rj^', which 
is obtained in an ensemble containing 20 random ini- 
tial orbit points {(a;i(0), ?/i(0), z,(0), 6*4(0)); i = 1,...,20} 
which are chosen with uniform probability in the range 
of a;,(0) e (-2,2), y,{Q) e (-16,0), z,(0) G (0,0.4), and 
6*^(0) e [0, 1). For the silent case of /dc = 0.21, F^^ grows 
up to the largest possible value of the derivative \dxi/d6\ 
along a trajectory and remains for all subsequent time. 
Thus, F^'' saturates for large N and hence the silent state 
has no phase sensitivity (i.e., it has smooth geometry). 
On the other hand, for the case of SN bursting, Py grows 
unboundedly with the same power 5, independently of i. 



r 



(0 

N 



(6) 



Here, the value of 5 ~ 1.74 is a quantitative character- 
istic of the phase sensitivity of the SN bursting attrac- 
tor for /dc = 0.24, and S is called the phase sensitivity 
exponent. For obtaining satisfactory statistics, we con- 
sider 20 ensembles for each /dc, each of which contains 



5 



20 randomly chosen initial points and choose the average 
value of the 20 phase sensitivity exponents obtained in 
the 20 ensembles. Figure ISjf) shows a plot of S versus 
A/dc (= ^dc — ^dc)- Note that the value of S monoton- 
ically increases from zero as I^c is increased away from 
the bursting transition point I^^ (~ 0.2236). As a result 
of this phase sensitivity, the SN bursting attractor has 
strange geometry leading to aperiodic complex bursting, 
as in the case of chaotic bursting attractor. 

Using the rational approximation to the quasiperiodic 
forcing (isl . , we explain the mechanism for the tran- 
sition from a silent to an SN bursting state. For the 
inverse golden mean, its rational approximants are given 
by the ratios of the Fibonacci numbers, u)k = Fk-i/Fk, 
where the sequence of {Fk} satisfies F^+i = Fk + Fk^i 
with Fq = and Fi = 1. Instead of the quasiperiodically 
forced system, we study an infinite sequence of periodi- 
cally forced systems with rational driving frequencies ujk- 
For each rational approximation of level fc, a periodically 
forced system has a periodic or a chaotic attractor that 
depends on the initial phase 6q of the external forcing. 
Then, the union of all attractors for different 6*0 gives the 
kth approximation to the attractor in the quasiperiodi- 
cally forced system. For this rational approximation of 
level fc, it is sufficient to change the initial phase h^ 
the interval [0,1/ Fk) in order to get all possible attract- 
ing sets, because the set of all values fills the whole 
interval [0, 1). 

We consider the rational approximation of level k = 7 
to the quasiperiodic forcing of A2 = 0.5. As shown in 
Fig. m^a) for I^c — 0.22, the rational approximation to 
a stable smooth torus (represented by a black curve), 
corresponding to a silent state, consists of stable orbits 
with period F7 (=13). Figure SKb) shows the magni- 
tude of the first stability multiplier Ai (with the largest 
magnitude) of the stable i^y-periodic orbits in the inter- 
val [0,1/^7). We note that |Ai| varies depending on 9, 
and all of its values are less than unity. Hence, all F7- 
periodic orbits for all 6 are stable. However, as I^c passes 
a threshold value /dc,7 (— 0.2209), the smooth torus be- 
comes broken and a nonsmooth bursting attracting set 
with Fi "gaps," where no stable orbits with period -FV 
exist, appears. An example is given in Fig. SJ^c) for 
^dc = 0.222. A magnified view of the 1st gap is given 
in the inset. For this case. Fig. |U[d) shows |Ai| of the 
F7-periodic orbits in the interval [0,1/^7). In the gap 
where 0.312 306 < F^O < 0.523 202, |Ai| is larger than 
unity, while in the remaining region of 9, |Ai| is less than 
unity. Thus, the ^V-periodic orbits in the gap become 
unstable via phase-dependent bifurcations (occurring at 
specific values of 9), and then chaotic bursting attractors 
fill the gap together with regular attractors with periods 
higher than Fj embedded in small windows. The bifur- 
cation diagram (i.e., plot of x vs. Fj9 in the i^V-times 
iterated Poincare map 

P^O for /dc = 0.222 IS given m 
Fig.Ufe). At both ends of the gap, the P7-periodic attrac- 
tor (denoted by a black solid curve) becomes unstable via 
a subcritical period-doubling bifurcation when it absorbs 
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FIG. 4: Investigation of transition from a smooth torus (cor- 
responding to a silent state) to an SN bursting attractor (cor- 
responding to an SN bursting state) in the rational approxi- 
mation of level k = 7 for Ai = 0.5, A2 = 0.5, and /i = 30 
Hz. In (a) and (c), projections of attractors onto the 6-x plane 
are given, (a) Stable smooth torus corresponding to the silent 
state for 7dc = 0.22. (b) Plot of |Ai| vs. FrO for Jdc = 0.22 [Ai: 
the first stability multiplier (with largest magnitude) of the 
fV-periodic orbit], (c) SN bursting attractor for Jdc ~ 0.222. 
A magnified view of the 1st gap is also given in the inset, (d) 
Plot of I All vs. FyO for /dc ~ 0.222. (e) Bifurcation diagram 
(i.e., plot of X vs. F76) for Jdc = 0.222 in the fV-times iter- 
ated Poincare map P^'' . An fV-periodic orbit (represented by 
a black solid line) becomes unstable via a subcritical period- 
doubling bifurcation, and then it is denoted by black solid 
circles, (f) Lyapunov-exponent diagram for Jdc = 0.222 in 
the Poincare map. For each randomly chosen initial point, 
we get the Lyapunov exponent by following the 10* Poincare 
maps after the transients of the 10'^ Poincare maps. 



an unstable orbit with doubled period 2F7 (represented 
by a gray short-dashed curve). Then, a jump to a chaotic 
bursting attractor [developed from the period-doubling 
cascade of the stable 2F7-periodic orbit (denoted by a 
gray solid curve)] occurs. Thus, in Fig. [UJc), the ra- 
tional approximation to the whole attractor consists of 
the union of the periodic component and the chaotic 
bursting component, where the latter occupies F^ gaps 
in 9. Figure [4jf) shows the Lyapunov-exponent diagram 
(i.e., plot of ai{9) vs. F79) for /dc = 0.222. (In the 
gap, chaotic bursting attractors with positive cti coexist 
along with periodic attractors with negative <ti embed- 
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ded in small windows.) The angle-averaged Lyapunov 
exponent (cti) [(•••) denotes the average over the whole 
6] is given by the sum of the "weighted" Lyapunov ex- 
ponents of the periodic and chaotic components, Ap and 
Ac, (i.e., (cTi) = Ap + Ac), where Ap(c) = Mp(c)(ai)p(c), 
and Mp(c) and (cri)p(c) are the Lebesgue measure in 
9 and the average Lyapunov exponent of the periodic 
(chaotic) component, respectively. Since the periodic 
component is dominant, the average Lyapunov exponent 
((cti) ~ —0.073) is negative. Hence, the rational approx- 
imation to the whole attractor in Fig.SJ^c) is nonchaotic. 
We note that Fig. UJ^c) resembles Fig. [3Kb), although the 
level fc = 7 is low. Increasing the level to k — 10, we 
confirm that the rational approximations to the whole 
attractor have Fk gaps (filled with chaotic bursting at- 
tractors) which appear via phase-dependent subcritical 
period-doubling bifurcations and their average Lyapunov 
exponents are negative. In this way, an SN bursting at- 
tractor appears in the case of quasiperiodic forcing, as 
shown in Fig. E^b). 

Finally, we characterize the bursting activity [alternat- 
ing between the silent phase and the active (bursting) 
phase] in both cases of SN and chaotic bursting states for 
A2 — 0.5. Figures [Sjal) and[5Ka2) show sequences of the 
interburst intervals [i.e., time interval between the first 
spikes in the neighboring bursts (i.e., active phases)] for 
the SN and chaotic bursting states when Idc = 0.24 and 
0.29, respectively. Both sequences are aperiodic complex 
ones. Histograms of the interburst intervals for these SN 
and chaotic bursting states are also given in Figs. E^bl) 
and[5Kb2), respectively. They are multimodal ones. As 
Idc is increased, heights of peaks for short interburst 
intervals increase, while those for longer ones decrease. 
Hence, as Idc is increased, the average interburst inter- 
val (IBI) ((••■) represents the average over a large num- 
ber of bursts) decreases {i.e., the mean bursting rate in- 
creases); (IBI) ~ 1029 ms and 549 ms for Idc = 0.24 
and 0.29, respectively. The active (bursting) phases are 
also characterized in terms of the bursting length (BL) 
(i.e., time interval between the first and last spikes in a 
burst) and the number of spikes (n) per burst. For our 
cases of the SN and the weakly chaotic burstings, some of 
the active phases are found to consist of only one spike 
(i.e., n=l), and hence their BLs are zero. Presence of 
these active phases with n=l implies "weak" firing ac- 
tivity in such active phases. (This is in contrast to the 
periodically forced case (of A2 = 0) where a direct tran- 
sition from a silent to a chaotic bursting state with all 
active phases composed of more than one spike occurs.) 
Sequences of both BL and n are shown in Figs. EKcl)- 
[5Kc2) and Figs. [5Kel)-[5Ke2). As in the case of the inter- 
burst intervals, they are aperiodic complex ones for both 
the SN and chaotic burstings. Likewise, their histograms 
are also multimodal, as shown in Figs. [5Kdl)-[SJd2) and 
Figs. [5Kfl)-[5Uf2). With increase in Idc, both the average 
bursting length (BL) and the average number of spikes 
(n) in a burst increase; (BL) ~ 23 ms and 39 ms, and 
(n) ~ 2.9 and 3.4 for Idc = 0.24 and 0.29, respectively. 
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FIG. 5: Characterization of SN {Idc = 0.24) and chaotic 
{Idc ~ 0.29) bursting states for Ai = 0.5, A2 = 0.5, and 
/i = 30 Hz. Sequences of interburst intervals (IBIs) for (al) 
SN and (a2) chaotic burstings; i represents the bursting index. 
Histograms of IBIs for (hi) SN and (b2) chaotic burstings. Se- 
quences of bursting length (BL) for (cl) SN and (c2) chaotic 
burstings. Histograms of BL for (dl) SN and (d2) chaotic 
burstings. Sequences of number of spikes (n) per burst for 
(el) SN and (e2) chaotic burstings. Histograms of n for (fl) 
SN and (f2) chaotic burstings. For histograms of IBI and 
BL, 200 equally-spaced bins are chosen in the range of IBI 
G (0, 2000) and BL e (0, 200). We get the number of bursts 
(A'^) in each bin among the total number of 5000 bursts (which 
are also used for the histogram of n). 



Thus, both the SN and chaotic bursting states exhibit 
aperiodic complex burstings, although their dynamics are 
different (one is chaotic and the other one is nonchaotic). 
We note that such aperiodic complexity results from the 
strange geometry of the SN and chaotic bursting states. 



III. SUMMARY 

We have investigated a dynamical transition from a 
silent state to a bursting state by varying the dc stimulus 
Idc in the quasiperiodically forced HR neuron. For this 
case of quasiperiodic forcing, a transition from a silent 
state to an SN bursting state (with negative Lyapunov 
exponent and positive phase sensitivity exponent) has 
been found to occur when Idc passes a threshold value. 
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With further increase in I^c, such an SN bursting state 
transforms to a chaotic bursting state (with a positive 
Lyapunov exponent). Thus, a new type of SN burst- 
ing states appear between the silent and chaotic bursting 
states as intermediate ones. This is in contrast to the 
periodically forced case where a dircict transition from 
a silent state to a chaotic bursting state occurs. Us- 
ing a rational approximation to the quasiperiodic forc- 
ing, the mechanism for the appearance of SN bursting 
states has been studied. Thus, a smooth torus, corre- 
sponding to a silent state, is found to transform to an SN 
bursting attractor through a phase-dependent subcritical 
period-doubling bifurcation. Both SN and chaotic burst- 
ing states have been characterized in terms of the inter- 
burst intervals, the bursting lengths, and the number of 
spikes per burst. As a result of their strange geometry. 



both bursting states are found to be aperiodic complex 
ones, although their dynamics are qualitatively different. 
Hence, we note that not only chaotic but also SN burst- 
ings may become dynamical origin of complex physiolog- 
ical rhythms which are central to life and ubiquitous in 
organisms. 
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